时间尺度重构EEMD-GRNN改进模型预测PM2.5的研究

合理构建PM_2.5浓度预测模型是科学、准确地预测PM_2.5浓度变化的关键。传统PM_2.5预测EEMD-GRNN模型具有较好的预测精度,但是存在过于关注研究数据本身而忽略其物理意义的不足。本研究基于南京市2014-2017年PM_2.5浓度时间序列数据,分析PM_2.5浓度多尺度变化特征及其对气象因子和大气污染因子的尺度响应,基于时间尺度重构进行EEMD-GRNN模型的改进与实证研究。南京市样本数据PM_2.5浓度变化表现为明显的天际尺度和月际尺度,从重构尺度（天际、月际）构建GRNN模型更具有现实意义;同时,PM_2.5对PM₁₀、NO₂、O₃、RH、MinT等因子存在多尺度响应效应,以其作为GRNN模型中的输入变量更具有时间序列上的解释意义。改进后的EEMD-GRNN模型具有更高的PM_2.5浓度预测精度,MAE、MAPE、RMSE和R²分别为6.17、18.41%、8.32和0.95,而传统EEMD-GRNN模型的模型有效性检验结果分别为8.37、27.56%、11.56、0.91。对于高浓度天（PM_2.5浓度大于100 μg/m³）的预测,改进模型更是全面优于传统EEMD-GRNN模型,MAPE为12.02%,相较于传统模型提高了9.03%。     

保持局域特征的多源海冰图像融合与分类

本文提出一种保持局域特征的多源海冰图像融合方法，并在此基础上进行海冰分类。本文提出的多源海冰图像融合方法包括保持空间局域融合和保持特征局域融合两方面。首先，通过学习得到投影矩阵和相似矩阵。投影矩阵将多源像素进行投影变换，得到保留像素空间局域特性的融合向量。相似矩阵度量像素特征间的相似性，通过拉普拉斯特征分解，得到保留像素特征局域相似性的融合向量。然后，将空间融合向量和特征融合向量进行像素综合，得到融合图像。在此基础上，本文设计一种滑动集成分类方法进行融合图像像素分类。提出的分类方法利用滑动集成的特点，在分类时增强刻画了海冰局域特性。由于本文的保持局域融合框架不仅刻画了海冰在物理空间中的邻接关系，而且考虑不同海冰类型的特征关系，因此其在多源图像（多光谱和合成孔径雷达（SAR）图像）的海冰分类任务中表现优异。实验结果表明本文提出的基于保持局域特征融合的多源海冰图像分类方法有效提升了海冰分类精度。     

Location and probability of shoal margin collapses in a sandy estuary

Channel bank failure, and collapses of shoal margins and beaches due to flow slides, have been recorded in Dutch estuaries for the past 200 years but have hardly been recognized elsewhere. Current predictions lack forecasting capabilities, because they were validated and calibrated for historic data of cross‐sections in specific systems, allowing local hindcast rather than location and probability forecasting. The objectives of this study were to investigate where on shoal margins collapses typically occur and what shoal margin collapse geometries and volumes are, such that we can predict their occurrence. We identified shoal margin collapses, generally completely submerged, from bathymetry data by analyzing digital elevation models of difference of the Western Scheldt for the period 1959–2015. We used the bathymetry data to determine the conditions for occurrence, specifically to obtain slope height and angle, and applied these variables in a shoal margin collapse predictor. We found 299 collapses along 300 km of shoal margin boundaries over 56 years, meaning that more than five collapses occur on average per year. The average shoal margin collapse body is well approximated by a 1/3 ellipsoid shape, covers on average an area of 34 000 m² and has an average volume of 100 000 m³. Shoal margin collapses occur mainly at locations where shoals take up a proportionally larger area than average in the cross‐section of the entire estuary, and occur most frequently where lateral shoal margin displacement is low. A receiver operating characteristic curve shows that the forecasting method predicts the shoal margin collapse location well. We conclude that the locations of the shoal margin collapses are well predicted by the variation in conditions of the relative slope height and angle within the Western Scheldt, and likely locations are at laterally relatively stable shoal margins. This provides hypotheses aiding the recognition of these features in sandy estuaries worldwide. Copyright © 2018 John Wiley & Sons, Ltd.     

Impacts of New Implementing Strategies for Surface and Model Physics Perturbations in TREPS on Forecasts of Landfalling Tropical Cyclones

To improve the ensemble prediction system of the tropical regional atmosphere model for the South China Sea (TREPS) in predicting landfalling tropical cyclones (TCs), the impacts of three new implementing strategies for surface and model physics perturbations in TREPS were evaluated for 19 TCs making landfall in China during 2014–16. For sea surface temperature (SST) perturbations, spatially uncorrelated random perturbations were replaced with spatially correlated ones. The multiplier f, which is used to form perturbed tendency in the Stochastically Perturbed Parameterization Tendency (SPPT) scheme, was inflated in regions with evident convective activity (f-inflated SPPT). Lastly, the Stochastically Perturbed Parameterization (SPP) scheme with 14 perturbed parameters selected from the planetary boundary layer, surface layer, microphysics, and cumulus convection parameterizations was added. Overall, all these methods improved forecasts more significantly for non-intensifying than intensifying TCs. Compared with f-inflated SPPT, the spatially correlated SST perturbations generally showed comparable performance but were more (less) skillful for intensifying (non-intensifying) TCs. The advantages of the spatially correlated SST perturbations and f-inflated SPPT were mainly present in the deterministic guidance for both TC track and wind and in the probabilistic guidance for reliability of wind. For intensifying TCs, adding SPP led to mixed impacts with significant improvements in probability-matched mean of modest winds and in probabilistic forecasts of rainfall; while for non-intensifying TCs, adding SPP frequently led to positive impacts on the deterministic guidance for track, intensity, strong winds, and moderate rainfall and on the probabilistic guidance for wind and discrimination of rainfall.     

Fengyun-4 Geostationary Satellite-Based Solar Energy Nowcasting System and Its Application in North China

Surface solar irradiance (SSI) nowcasting (0–3 h) is an effective way to overcome the intermittency of solar energy and to ensure the safe operation of grid-connected solar power plants. In this study, an SSI estimate and nowcasting system was established using the near-infrared channel of Fengyun-4A (FY-4A) geostationary satellite. The system is composed of two key components: The first is a hybrid SSI estimation method combining a physical clear-sky model and an empirical cloudy-sky model. The second component is the SSI nowcasting model, the core of which is the derivation of the cloud motion vector (CMV) using the block-matching method. The goal of simultaneous estimation and nowcasting of global horizontal irradiance (GHI) and direct normal irradiance (DNI) is fulfilled. The system was evaluated under different sky conditions using SSI measurements at Xianghe, a radiation station in the North China Plain. The results show that the accuracy of GHI estimation is higher than that of DNI estimation, with a normalized root-mean-square error (nRMSE) of 22.4% relative to 45.4%. The nRMSE of forecasting GHI and DNI at 30–180 min ahead varied within 25.1%–30.8% and 48.1%–53.4%, respectively. The discrepancy of SSI estimation depends on cloud occurrence frequency and shows a seasonal pattern, being lower in spring–summer and higher in autumn–winter. The FY-4A has great potential in supporting SSI nowcasting, which promotes the development of photovoltaic energy and the reduction of carbon emissions in China. The system can be improved further if calibration of the empirical method is improved.     

Comparison of different land-surface perturbation methods in short-range ensemble forecasts

为了比较不同陆面扰动方法对短时集合预报的影响,本研究设计了陆面模式扰动实验(LSMPE),初始土壤湿度扰动实验(ISMPE),陆面-大气耦合系数扰动实验(LCCPE)以及大气扰动对照实验(GEFSPE).结果表明,在三组陆面扰动实验中,LSMPE能代表最大的不确定性且误差最小;ISMPE的离散度要比LCCPE稍大,但是耦合系数的减弱能有效降低区域平均降水量.然而,所有陆面扰动实验产生的离散度均小于GEFSPE,这表明未来有必要将陆而扰动与大气扰动相结合,以使得模式能代表更多的不确定性.     

基于CMIP6模式优化集合平均预估21世纪全球陆地生态系统总初级生产力变化

利用国际耦合模式比较计划第六阶段(CMIP6)中18个地球系统模式总初级生产力(GPP)模拟数据,基于传统的多模式集合平均(MME)和可靠集合平均方法(REA),在4个未来情景(SSP1-2.6,SSP2-4.5,SSP3-7.0和SSP5-8.5)下预估了21世纪全球陆地生态系统GPP的变化量,并分析了GPP变化的驱动因子。研究结果表明：在4个未来情景下,基于REA方法预估的全球陆地生态系统年GPP在未来时期(2068—2100年)比历史时期(1982—2014年)分别增长了(14.85±3.32)、(28.43±4.97)、(37.66±7.61)和(45.89±9.21)Pg C,其增量大小和不确定性都明显低于MME方法。在4个情景下,大气CO₂浓度增长对GPP变化的贡献最大,基于REA方法计算的贡献占比分别为140%、137%、115%和75%;除SSP5-8.5(24%)外,其他情景下升温均导致全球陆地生态系统GPP降低(-42%、-37%、-16%),部分抵消了CO₂施肥效应的正面贡献。温度的影响存在纬度差异：升温在低纬度地区对GPP有负向贡献,在中高纬度地区为正向贡献。降水和辐射变化对GPP变化的贡献相对较小。     

概率密度分布法在“一带一路”地区极端月气温评估和预估中的应用

“一带一路”地区人口众多,气候类型复杂,亟待加强区域气候变化风险的认识。文中将该区分成10个区域,基于第五次耦合模式比较计划（CMIP5）中的31个全球模式模拟结果,应用概率密度分布（PDF）方法评估历史阶段（1986—2005年）各模式模拟暖月和冷月气温的能力,挑选并建立较优模式集合,用以预估21世纪中叶（2041—2060年）和21世纪末（2081—2100年）的极端月气温。结果表明,模式对观测中冷月气温距平PDF的模拟水平整体较暖月高。与多模式平均以及中位值相比,较优模式集合方法更适于极端暖/冷月气温的评估。在中等排放RCP4.5情景下,与低纬度地区相比,较优模式模拟中高纬地区未来极端暖/冷月气温的增温幅度的不确定性范围较大。21世纪中叶和21世纪末较优模式模拟的极端暖月气温在地中海增幅整体最大,东南亚增幅整体最小。对较优模式集合预估的极端冷月气温而言,无论是21世纪中叶还是世纪末,北欧增幅整体最大,东南亚增幅整体最小。     

基于卡尔曼滤波的中国区域气温和降水的多模式集成预报

利用TIGGE资料集下欧洲中期天气预报中心(ECMWF)、日本气象厅(JMA)、美国国家环境预报中心(NCEP)、中国气象局(CMA)和英国气象局(UKMO)5个模式预报的结果,对基于卡尔曼滤波的气温和降水的多模式集成预报进行研究。结果表明,卡尔曼滤波方法的预报效果优于消除偏差集合平均(BREM)和单模式的预报,但是对于地面气温和降水,其预报效果也存在一定的差异。在中国区域2 m气温的预报中,卡尔曼滤波的预报结果最优。而对于24 h累积降水预报,尽管卡尔曼滤波在所有量级下的TS评分均优于BREM,但随着预报时效增加,其在大雨及以上量级的TS评分跟最佳单模式UKMO预报相当,改进效果不明显。卡尔曼滤波在地面气温和24 h累积降水每个预报时效下的均方根误差均最优,预报效果更佳且稳定。     

A CRPS-Based Spatial Technique for the Verification of Ensemble Precipitation Forecasts

Traditional precipitation skill scores are affected by the well-known“double penalty”problem caused by the slight spatial or temporal mismatches between forecasts and observations. The fuzzy (neighborhood) method has been proposed for deterministic simulations and shown some ability to solve this problem. The increasing resolution of ensemble forecasts of precipitation means that they now have similar problems as deterministic forecasts. We developed an ensemble precipitation verification skill score, i.e., the Spatial Continuous Ranked Probability Score (SCRPS), and used it to extend spatial verification from deterministic into ensemble forecasts. The SCRPS is a spatial technique based on the Continuous Ranked Probability Score (CRPS) and the fuzzy method. A fast binomial random variation generator was used to obtain random indexes based on the climatological mean observed frequency, which were then used in the reference score to calculate the skill score of the SCRPS. The verification results obtained using daily forecast products from the ECMWF ensemble forecasts and quantitative precipitation estimation products from the OPERA datasets during June-August 2018 shows that the spatial score is not affected by the number of ensemble forecast members and that a consistent assessment can be obtained. The score can reflect the performance of ensemble forecasts in modeling precipitation and thus can be widely used.